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Abstract. We present a parametrization of the observed enhancement in the transverse electron quasielas- 
tic (QE) response function for nucleons bound in carbon as a function of the square of the four momentum 
transfer (Q 2 ) in terms of a correction to the magnetic form factors of bound nucleons. The parametriza- 
tion should also be applicable to the transverse cross section in neutrino scattering. If the transverse 
enhancement originates from meson exchange currents (MEC), then it is theoretically expected that any 
enhancement in the longitudinal or axial contributions is small. We present the predictions of the "Trans- 
verse Enhancement" model (which is based on electron scattering data only) for the v^, differential and 
total QE cross sections for nucleons bound in carbon. The Q 2 dependence of the transverse enhancement 
is observed to resolve much of the long standing discrepancy in the QE total cross sections and differential 
distributions between low energy and high energy neutrino experiments on nuclear targets. 

PACS. 13.15.+g Neutrino interactions - 25.30.Pt Neutrino scattering - 13.40.Gp Electromagnetic form 
factors 



1 Introduction 

A reliable description of the neutrino {v^) and antineu- 
trino (v^j) quasielastic (QE) and inelastic scattering pro- 
cesses (particularly on nuclear targets) is essential for pre- 
cision studies of oscillation [TJ[2] parameters such as 
mass splitting and mixing angles. In addition to model- 
ing the v^v^ cross sections [3J, a reliable model of the 
hadronic final states is needed because the hadronic en- 
ergy response of detectors is not the same for protons, 
neutrons, pions, photons, and nuclear fragments. Prescrip- 
tions which can be readily incorporated into existing i/ M 
Monte Carlo generators [1] are preferable. 

Models which assume that QE scattering on nuclear 
targets can be described in terms of scattering from inde- 
pendent nucleons bound in a nuclear potential (e.g. Fermi 
gas [5 or spectral functions) do not provide an adequate 
representation of measured differential and total QE cross 
sections for low energy (ss lGeV) Vn scattering on nucle- 
ons bound in carbon[Sl|7] (MiniBooNE) and oxygen[Sl[H] 
(K2K and T2K). The measured QE total cross sections 
are 20% larger than the model and the differential distri- 
butions in Q 2 are also inconsistent. The vector and axial 
form factors that are used in independent nucleon models 
are the free nucleon form factors extracted from electron 
and Vfi, scattering data on hydrogen and deuteriumjlO . 

Although there are more sophisticated calculations of 
quasielastic scattering (e.g. relativistic distorted-wave im- 
pulse approximationpT)). it is the simple independent nu- 



cleon model that has been implemented in the currently 
available neutrino cross section Monte Carlos 0]. 

This disagreement between the measured low energy 
differential and total QE cross sections on nuclear tar- 
gets and the predictions from the independent nucleon 
model has been attributed to an incomplete description of 
nuclear effects. These additional nuclear effects have been 
parametrized as an ad-hoc change in the the axial form 
factor mass parameter from the value measured for free 
nucleonsQU] (Al{ ree = 1.014 ± 0.014 GeV) to M e J f = 

1.20 ± 0.12 GeV (K2K) and M e J s = 1.23 ± 0.20 GeV 
(MiniBooNE). 

A recent analysis [7] of newly published differential QE 
cross sections from MiniBooNE (for nucleons bound in 
carbon) yields larger values of = 1.350 ±0.066 GeV 

in the Fermi gas model and M e J f = 1.343 ± 0.060 GeV 
in the spectral function model. In that analysis the free 
nucleon value M^ ee — 1.014 GeV is excluded at the con- 
fidence level greater than 5a (standard deviations). 

Figure II] shows the world's data [10] for the nucleon 
axial form factor (Fa(Q 2 )) extracted from QE v^, scat- 
tering on hydrogen and deuterium. Here, the data for 
Fa{Q 2 ) are shown as a ratio to a nominal dipole Fa(Q 2 ) = 
G^(Q 2 )= -^ga with M A = 1-015 GeV. On the left 

side we show the values extracted from ia, , experiments 
on hydrogen and deuterium and on the right side we show 
the values extracted from pion electro-production data on 
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Fig. 1. (a) F A (Q 2 ) extracted from ^-deuterium data[lO] divided by Go(Q 2 ) with M A = 1.015 GeV. (b) F A (Q 2 ) from 
pion electroproduction (corrected for hadronic effects) divided by Gp(Q 2 ) with Ma = 1.015 GeV. Thin solid line - duality 
(10] : Short-dashed line - Fa{Q 2 )a2=v2- Dashed-dot line - constituent quark model; Thick solid red line 
i_i_ n2'/» J2 ^2 with Ma = 1-35 GeV. The horizontal scale on top is Q 2 . The horizontal scale on the bottom 

->2 , 



based fit from reference 
Fa(Q 2 ) = G£(Q 2 )=- ( 

is the target mass scaling variable £ for elastic scattering (x — 1) 
nucleon mass. 



Here £ = 



r = Q J /4M- i , and M is the average 



hydrogen. The average of the measurements of Ma from 
z/ p , experiments on hydrogen and deuterium of Ma'"" 
= 1.016±0.026 GeV is in agreement with the average value 
of Mj^ on —1.014:± 0.016 GeV extracted from pion electro- 
production experiments on hydrogen (after corrections for 
hadronic effects). The average of the v^Vn an d electro- 
production values is [TU] M% orld - av = 1.014 ±0.014 GeV. 
The thin solid line is a duality based parametrization|10| 
of possible deviations from the dipole form. The dashed- 
dot line is the prediction of a constituent quark model [II] 
and the short-dashed line is the expectation for Fa(Q 2 ) if 
the vector and axial-vector structure functions are equal 

l } a jQelastic— vector -» \ <Qelastic— axial\ 

(eg. vv 2 — vv 2 )■ 

It is clearly observed that a dipole axial form factor 
with Ma = 1-35 GeV (thick solid red line) is inconsistent 
with the measurements on hydrogen and deuterium. 

It has been assumed that an " effective" axial mass 
provides an adequate description of the missing nuclear 
corrections. However, a large increase in the axial form 
factor of bound nucleons is contrary to theoretical expec- 
tations that Ma in nuclear targets should be smaller 15 
than (or the same|16j) as in deuterium. 

Additionally, the low energy neutrino data appear to 
be in disagreement with higher energy neutrino experi- 
ments on nuclear targets. At high neutrino energies, the 
total and differential QE cross sections on nuclear targets 
are consistent with models which assume that the scat- 
tering is on independent nucleons with free nucleon form 
factors. For example, Ma of 0.979 ± 0.016 GeV has been 
extracted from a global analysis [17] of the differential dis- 



tributions and total QE cross sections measured in all high 
energy experiments on nuclear targets. 

Recent measurements of the differential and total QE 
cross section for nucleons bound in carbon by the NOMAD [PE] 
collaboration for i/^, energies above 4 GeV are also con- 
sistent with models which assume that the scattering is 
from independent nucleons with free nucleon form factors. 
The NOMAD analysis yields a value of M A (1.05 ± 0.02 ± 
0.06 GeV) 

Therefore, the results of the higher energy and low 
energy experiments on nuclear targets appear to be 
inconsistent with each other. 

In this communication we investigate the transverse 
enhancement observed in QE electron scattering experi- 
ments on nuclear targets. We obtain a parametrization of 
the enhancement and investigate its implication for v^, 
scattering. We show that the Q 2 dependence of the trans- 
verse enhancement resolves much of the discrepancy be- 
tween the low energy and high energy neutrino experi- 
ments, in addition to obviating the need for an ad-hoc 
nuclear modification to Ma- 



2 Electron-nucleon scattering 

The differential cross section for scattering of an unpo- 
larized charged lepton with an incident energy Eq, final 
energy E' and scattering angle 9 can be written in terms 
of the structure functions T\ and T 2 as: 

^ {Eo ,E>,e) = ^cos 2 (e/2) 
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x [? 2 (x, Q 2 )/u + 2 tan 2 (0/2)71 (z, Q 2 )/M] 

where a is the fine structure constant, M is the nucleon 
mass, v = Eq — E' is the energy of the virtual photon 
which mediates the interaction, Q 2 — AEqE' sin 2 (6'/2) is 
the invariant four-momentum transfer squared, and x = 
Q 2 /2Mv is the Bjorken scaling variable. We define T 2 = 
vW 2 , T\ = MWi (and for scattering IF3 — vW$). 

Alternatively, one could view this scattering process 
as virtual photon absorption. Unlike the real photon, the 
virtual photon can have two modes of polarization. In 
terms of the cross section for the absorption of transverse 
(ox) and longitudinal (ctl) virtual photons, the differen- 
tial cross section can be written as, 



d 2 o 
dQdE 



7 = r[a T {x,Q 2 ) + ec7 L {x,Q 2 )] 



where, 



r 



K 



aKE' 
4w 2 Q 2 E 

1 + 2(1 + 

2Mv - Q 2 
2M 



1 - e 
Q 2 



-U+\ C 2 )tan2 2 



(1) 

(2) 

(3) 
(4) 



The quantities r and e represent the flux and the de- 
gree of longitudinal polarization of the virtual photons 
respectively, which the quantity R is defined as the ratio 
cl/cT) a n d is related to the structure functions by 



R{x,Q 2 ) = 



a 'I 



T 2 

2xT\ 



(1 



4M 2 z 2 . 



- 1 = 



2xTi 



(5) 



where JF L is called the longitudinal structure function. The 
structure functions are expressed in terms of 07, and ot 
as follows: 



T-2 



MK 
Air z a 
vK{ol + or) 



^l(x,Q 2 



47t 2 q;(1 
T 2 ( 1 



r 4M 2 x 2 

AM 2 x 2 



2xF x 



or. 



IxTx = T 2 (l + ) - Fl(x, Q 2 ). 

In addition, 2xT\ is given by 



(6) 
(7) 

(8) 
(9) 



2xT 1 (x,Q 2 ) = T 2 (x,Q' z ) 



2 , 1 + AM 2 x 2 /Q 2 



or equivalently 



Wi(x,Q 2 ) = W 2 {x,Q 2 



l + R(x,Q 2 ) 



1 + v 2 IQ 2 
l + R{x,Q 2 ) 



In the case of elastic scattering from free nucleons (x = 
Q 2 /2Mv=l) the structure functions are related to the nu- 
cleon form factors by the following expressions [19 : 



yyelastic 



Q 2 



2\|2 
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w. 
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In 
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2p 
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dasii 



. g( _ Q 2 , [G Ep {Q 2 )} 2 + t[G Mp (Q 2 )} 2 
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§( _ Q 2 AG En (Q 2 )] 2 +r[G Mn {Q 2 )] 2 
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R 



elastic 



(x = l,Q 2 ) 



-elastic 



T elastic 



4M 2 ( G\ 



Q 2 \G 2 m 



Here, r = Q 2 /AM 2 n , where M p .„ are the masses of proton 
and neutron. Therefore, Gmp and Gmu contribute to the 
transverse virtual photo-absorption cross section, and G e p 
and Geu contribute to the longitudinal cross section. 



3 Nucleon form factors 

The nucleon electromagnetic form factors are best de- 
scribed by the BBBA2{)07 2 5 duality based parametrization 10J 
The deviations from the dipole form factors are parametrized 
by multiplicatives functions Ajv(£) for each of the proton 
and neutron form factors (A Ep ((, p ), A Mp (£ p ), A En (£ n ), 
and Am™ (£"))■ Here, -A/v(£) = 1 for pure dipole form fac- 
tors. The variable £ is the target mass scaling variable for 
clastic scattering (x — 1), where 



(1 + y/l + l/r p , n ) ' 



and T p>n = Q 2 /AM 2 n . Here M Pi „ are the proton (0.9383 
GeV/c 2 ) and neutron (0.9396 GeV/c 2 ) masses, respec- 
tively. 

g v d (Q 2 ) 



(1 + Q 2 /M 2 Y 

Ge P (Q 2 ) = AEp~dipole{£ P ) X GYf^Q 2 ) 

G En (Q 2 ) = Af n (C) x G Ep {Q 2 ) x (r^T 

Gm P (Q 2 )/H P — AM P -di P ole{i P ) X GY)(Q 2 ) 

G Mn (Q 2 )/u n = Af In (C) x G Mp (Q 2 )/u p 

Here fi p = 2.7928, u n = -1.913, and M$ = 0.71 GeV 2 
(My = 0.8426 GeV). The parameters for the multiplica- 
tive functions An(£) which describes the ratio to dipole 
are given in reference [TU]. The parametrizations are com- 
pared to experimental data in Figure [2] 
For the axial form factor we use 



F A (Q 2 ) = Gi(Q 2 ' 



flu 



(1 + Q 2 /M 2 ) 2 
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Fig. 2. Ratios of G Ep (a), G Mp /^ p (b), G En (c) and G Mn /fJ.n (d) to G V D = (1+Q2 ) M 2 p with My = 0.8426 GeV. The short- 
dashed line in each plot is the old Kelly [12] parameterizations (old Galster[l3] for G E n)- The solid line is the BBBA0725 and 
the long-dashed line is BBBAQT^ parametrizations[TU], respectively. The values of £ and the corresponding values of Q 2 are 
shown on the bottom and top axis. 



where g A = -1.267, and M A = 1-014 ± 0.014 GeV is the 
axial mass for free nucleons. 

The ratio of longitudinal and transverse cross sections 
for free nucleons is given by: 



R 



elastic _ 4MVj4 A Ep-dipole _ 0.481 A Ep-dipole 

r> 2 a 2 n 2 a 2 

^Mp-dipole ^Mp-dipole 



R 



' i ■ t * I ! clastic 



In the dipole approximation with Geu — 



R 



ilasti' 
deutei 



4,M 2 /(fi 2 p + nl) 0.328 



Q 2 



Q 2 



(10) 



4 Quasielastic electron scattering from 
nuclear targets 

For electron-nucleon and muon-nucleon scattering, scat- 
tering from free nucleons (with no pions in the final state) 



is called elastic scattering, and scattering from nucleons 
bound in a nuclear target (with no pions in the final 
state) is called QE scattering because the scattering is 
from quasi-free nucleons. 

For charged-current -nucleon and ^-nucleon scat- 
tering (with no pions in the final state), the term QE scat- 
tering is used to describe scattering from either free or 
bound nucleons because the neutrino is transformed to a 
final state muon. For neutrino processes, the term elastic 
scattering is only used when there is a neutrino in the final 
state. 

Studies of QE electron scattering on nuclear targets in- 
dicate that only the longitudinal part of the QE cross sec- 
tion can be described in terms of a universal response func- 
tion of independent nucleons bound in a nuclear potential|20j 
(and free nucleon form factors). In contrast, a significant 
additional enhancement with respect to the model is ob- 
served in the transverse part of the QE cross section. 

The enhancement in the transverse QE cross section 
has been attributed to meson exchange currents (MEC) in 
a nucleus [2Tjll2"Tll2"2"ir23U24) . Meson exchange currents orig- 
inate from nucleon-nucleon correlations (predominantly 
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neutron- proton). The final state for the MEC process 
can include one or two nucleons. If no final state pions are 
produced, the process is considered as an enhancement of 
the QE cross section. If one or more final state pions are 
produced, the process enhances the inelastic cross section. 

Within models of meson exchange currents the en- 
hancement is primarily in the transverse part of the QE 
cross section, while the enhancement in the longitudinal 
QE cross section is small (in agreement with the electron 
scattering experimental data). The conserved vector cur- 
rent hypothesis (CVC) implies that the corresponding vec- 
tor structure function for the QE cross section in , 
scattering can be expressed in terms of the structure func- 
tions measured in electron scattering on nuclear targets. 
Therefore, there should also be a transverse enhancement 
in neutrino scattering. 

In addition, for some models of meson exchange currents 
the enhancement in the axial part of v^, QE cross sec- 
tion on nuclear targets is also small. Therefore, the axial 
form factor for bound nucleons is expected to be the same 
as the axial form factor for free nucleons. 



4.1 Measuring the transverse enhancement at low Q 2 

The longitudinal response scaling functions extracted by 
Donnely et. al.[20] for different momentum scales and dif- 
ferent nuclei (A=12 ,40 and 56) are essentially described 
by one universal curve [20 which is a function of the nu- 
clear scaling variable ip' only. The function peaks at ip'=0 
and ranges from ip' = —1.2 to ip' = 2. In contrast, the 
transverse response scaling function is larger and increases 
with momentum transfer. The response function of the 
transverse enhancement excess is shifted to higher ip' and 
peaks at ip' « 0.2. 

Carlson et. alJJS] uses the measured longitudinal and 
transverse response functions to extract the ratio (TZt) of 
the integrated response functions for the transverse and 
transverse components of the QE response functions for 
values of ip' < 0.5 and ip' < 1.2. 

For nucleons bound in carbon, the ratios for ip' < 0.5 
are 1.2, 1.5, 1.65 for values of the 3-momcntum transfer 



to the ratio for ip' < 0.5 by the following expression: 

K T (all - ip') = 1 + 1-18 [n T (ip' < 0.5) - 1] 

We obtain TZ T (all-ip') values of 1.24±0.1, 1.59±0.1, and 
1.77 ± 0.1 for Q 2 values of 0.09, 0.15, and 0.33 (GeV/c) 2 , 
respectively. We use the difference in the measured values 
of TZt for ip' < 0.5 and ip' < 1.2 as an estimate of the 
systematic error. Since the longitudinal response function 
is equal to the response function for independent nucleons, 
the ratio TZt(o,11 — ip') is equivalent to the ratio of the 
integrated transverse response function in a nucleus to the 
response function for independent nucleons (as a function 
ofQ 2 ). 

The values of TZt extracted from the data of from Carl- 
son et al are shown as a function of Q 2 (black points) in 
g|f gure |] 

Transverse Enhancement Carbon 12 



Q 



qs of 0.3, 0.5, and 0.6 GeV/c, respectively (q 2 
where v = Q 2 /2M at the QE peak). 

The ratios for ip' < 1.2 are 1.25, 1.6, 1.8 for q^ values of 
0.3, 0.5, and 0.6 GeV, respectively. (These correspond to 
Q 2 values of 0.09, 0.15, and 0.33). At higher values of ip' 
the transverse response functions include both QE scat- 
tering and pion production processes (e.g. A production 
with Fermi motion). 

Therefore, we use the measured values of TZt for ip' < 
0.5, where the contribution from pion production process 
is small, and apply correction to extract the ratio for the 
entire range of ip' , as described below. 

The excess transverse response function peaks at ip' s» 
0.2, while the longitudinal response function peaks at ip' = 
0. A fit of an asymmetric gaussian to the longitudinal 
response function indicates that the TZt values for the 
total response functions integrated over all ip' are related 
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Fig. 3. The transverse enhancement ratio (TZt) as a func- 
tion of Q 2 . Here, TZt is ratio of the integrated transverse re- 
sponse function for QE electron scattering on nucleons bound 
in carbon divided by the integrated response function for in- 
dependent nucleons. The black points are extracted from Carl- 
son et <>/|23|. and the blue bands are extracted from a fit|29| 
to QE data from the JUPITER 25] experiment (Jlab exper- 
iment E04-001). The curve is a fit to the data of the form 
TZ T = 1 + AQ 2 e~ Q /B . The dashed lines are the upper and 
lower error bands. 



4.2 Measuring the transverse enhancement at high Q 2 

The technique of using the ratio of longitudinal and trans- 
verse QE structure functions to determine the transverse 
enhancement in the response functions for QE scattering 
is less reliable for Q 2 > 0.5 (GeV/c) 2 , because at high 
values of Q 2 the longitudinal contribution to the QE cross 
section is small (as illustrated in equation 10). 

Since the transverse cross section dominates at large 
Q 2 one can extract the transverse enhancement by com- 
paring the measured QE cross sections to the predictions 



6 



A. Bodek, H.S. Budd and M. E. Christy: Neutrino Quasielastic Scattering on Nuclear Targets 



Preliminary E04-001 , E = 1 .204-, - 28.01 1 



Preliminary E04-001 , E = 1 .204, = 70.01 1 



I 7 000 

1 oooo 
sooo 

6000 
4000 




0.6 O 



W z (GeV z ) 

Preliminary E04-00 1 , E = 4.62 9, = 10.661 



Q 2 = 0.68 (GeV/c) 2 
e = 0.98 
R T = 1.7 



^1 



Total 
QE 

Inelastic 
TE 

QE transverse 
QE Longitudinal 




*2 300 



■E 200 



Q 2 = 1 


.0 (GeV/c) 2 


e — O. 


44 


R T = 1 


.35 




0.5 O 



Preliminary E04-001 , E = 2.347, = 30.01 1 




W'(GeV ) 



Fig. 4. Samples of fits[29| to preliminary electron scattering 
data from the JUPITER collaboration[25 (Jefferson Lab ex- 
periment E04-001) on a carbon target. Shown are the contri- 
butions from the transverse QE (solid pink), longitudinal QE 
(dashed pink), total QE (solid red), inelastic (pion production) 
processes (solid green), and a transverse excess (TE) contribu- 
tion (solid black line). Top: Q 2 = 0.3 GeV/c 2 at the QE peak. 
Bottom: Q 2 = 0.68 GeV/c 2 at the QE peak. 

of the independent nucleon model directly. However, be- 
cause there is overlap between pion production processes 
and QE scattering, the contribution from pion production 
processes should be accounted for in the extraction pro- 
cess. 

We extract the transverse enhancement at higher val- 
ues of Q 2 from a fit to both existing electron scattering 
data on nuclei and preliminary data from the JUPITER 
collaboration 25J (Jefferson lab experiment E04-001). The 
fit (developed by P. Bosted and V. Mamyan) provides a 
description of inclusive electron scattering cross sections 
on a range of nuclei with A > 2. It is an extension of fits to 
the free proton [57] and deuteron [35] and was utilized for 
calculations of the radiative corrections for the JUPITER 
analysis [23]. Experiment E04-001 was designed to provide 
separations of the longitudinal and transverse structure 
functions from a range of nuclei. These data, therefore, 
provides a significant constraint on this separation in both 
the quasi-elastic and resonance regions, which is of criti- 



Fig. 5. Same as Fig. |4j Top: Q 2 = 1.0 GeV/c 2 at the QE peak. 
Bottom: Q 2 = 1.1 GeV/c 2 at the QE peak). 

cal importance for the current study. A brief description 
of the fit is given in [53] , which also provides plots of the 
fit residuals to the data sets utilized. 

The inclusive fit is a sum of four components: 

— The longitudinal QE contribution calculated for inde- 
pendent nucleons (smeared by Fermi motion in car- 
bon) 

— The transverse QE contribution calculated for inde- 
pendent nucleons (smeared by Fermi motion in car- 
bon) 

— A transverse excess (TE) contribution 

— The contribution of inelastic pion production processes 
(smeared by Fermi motion in carbon). 

The QE model used in the Bosted-Mamyan fit is the 
super-scaling model 26 of Sick, Donnelly, and Maieron. 

Figures [3] and [5] show samples of Bosted-Mamyan fits 
to preliminary electron scattering data from JUPITER 
on a carbon target. Shown are the contributions from 
the transverse QE (solid pink), longitudinal QE (dashed 
pink), total QE (solid red), inelastic pion production pro- 
cesses (solid green), and a transverse excess (TE) contri- 
bution (solid black line). 

We extract the transverse enhancement ratio as a func- 
tion of Q 2 by integrating the various contributions to the 
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fit up to W 2 = 1.5 GeV 2 . Here 

Q -^transverse H~ TE 



TZ T 



QE, 



ransverse 



We assign a conservative systematic error to TZt to ac- 
count for the possibility that a fraction of the transverse 
excess events may be produced with a pion in the final 
state. 

Figure [3] shows the values of TZt as a function of Q 2 . 
The black points are extracted from Carlson et al[2'6\. 
and the higher Q 2 blue bands are from the fit to QE 
data from the JUPITER collaboration^] . The data are 
parametrized by the expression: 

TZ T = 1 + AQ 2 e- Q ' 2 ' B 

with A = 6.0 and B = 0.34 (GeV/c) 2 . The electron scat- 
tering data indicates that the transverse enhancement is 
maximal near Q 2 = 0.3 (GeV jc) 2 and is small for Q 2 
greater than 1.5 (GeV/c) 2 . The upper error band is given 
by A = 6.7 and B = 0.35 (GeV/c) 2 , and the lower error 
band is given by A = 5.3 and B = 0.33 (GeV/c) 2 . This 
parametrization is valid for carbon (A=12) and higher A 
nuclei. 



Fermi Gas Correction 



o 1 



o 0.8 



y 0.6 



0.4 



0.2 



















































+ — Ferrhi Gas Correction 






iii, 


, , , , 











O 0.05 O.I O.I 5 0.2 0.25 0.3 

Q 2 (GeV/c) 2 



5.1 The "Independent Nucleon (Ma=1.014)" baseline 
model 

In modeling , QE scattering on nuclear targets we use 
BBBA200725 free nucleon electromagnetic form factors 
(with My = 0.71), and a dipole axial form factor with 
¥4 = 1.014 GeV. We apply Pauli blocking corrections 
to the differential QE cross section, as implemented in the 
NEUGEN Monte Carlo ^ . The Pauli blocking factor as a 
function of Q 2 is shown in figure [6] We do not apply Fermi 
motion corrections since we only study the total integrated 
QE cross section. We refer to this baseline model, which is 
shown as orange dotted lines on plots, as the "Independent 
Nucleon (Af A =1.014)" model. 



5.2 The "Transverse Enhancement" model 

We use our parametrization of TZt to modify Gmp and 
Gmu for bound nucleons as follows. First, we assume that 
the enhancement in the transverse QE cross section modi- 
fies Q\ l — Gmp — Gmu for nucleons bound in carbon with 
a form given by 



G 



nuclear 



Alp 
Qnucl 



nuclear 
Mn 



(Q 2 ) = G Mp (Q 2 ) x Vl + AQ 2 e-Q 2 /z 
(Q 2 ) = G Mn (Q 2 ) x Vl + AQ 2 e~Q 2 / B . 



In all of the studies we keep Ge p (Q 2 ), Geu(Q 2 ) and Fa(Q 2 ) 
for bound nucleons the same as for free nucleons. The 
transverse enhancement leads to an enhancement in the 
structure functions W® elasUc , w Q^u c and ^Qelastic 

The expressions for the f^,^ differential QE cross sec- 
tions are given in the Appendix. We also apply Pauli 
blocking as a function of Q 2 as shown in figure 6 We re- 
fer to this model as the "Transverse Enhancement" model. 
The predictions on the plots for the " Transverse Enhance- 
ment Model" are shown with solid red lines. The error 
bands are shown as dotted dashed red lines. The ratio 
of calculated quantities for the " Transverse Enhancement 
model" divided by "Independent Nucleon (ma=1.014) are 
also shown as solid red lines. 



Fig. 6. The Fermi suppression factor (Pauli blocking) used 
in our studies as a function of Q 2 . We use the Pauli blocking 
factor which is implemented in the NUEGEN Monte Carlo[l]. 



5.3 The "Larger M A (M A =1.3)" model 



5 Consequences for 
scattering on carbon 



charged-current QE 



We assume that there is a corresponding transverse en- 
hancement in the i/u, v^. QE cross sections on nuclear tar- 
gets. Although motivated by MEC, the analysis is model 
independent since the parameters are taken from electron 
scattering data. 

In the rest of this paper, the terms cross sections and 
differential distributions refer to scattering from nucleons 
bound in carbon. 



ons with 



THnuclear / 



Since low energy neutrino experiments have used an ad- 
hoc w 1.3 GeV to account for additional nuclear 
effects, we also compare our results to the differential and 
total QE cross sections calculated for independent nucle- 
1.3 GeV in the following expression: 

(q2) = ( i + cW 

For this model, we use the electromagnetic form factors 
for free nucleons, and apply Pauli blocking as described 
above. We refer to this model, which is shown as dashed 
blue lines on plots, as the "Larger Ma (Ma=1.3)" model. 
The ratio of calculated quantities for the "Larger Ma 
(Mj4=1.3)" model divided by the predictions of the "In- 
dependent Nucleon (71^4=1.014) model are also shown as 
dashed blue lines. 
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Fig. 7. The QE differential cross section (da/dQ 2 ) as a tunc- Fig. 8. Same as figure 7 for v^,Vn energies of 3.0 GeV (max- 



tion of Q for v^, energies of 1.0 GeV (maximum accessible imum accessible Q 
Qmax = 1-3 (GeV/c) 2 ). Here, the orange dotted line is the 
prediction of the "Independent Nucleon (Ma=1.014)" model. 
The blue dashed line is the prediction of the the "Larger Ma 
(A/a = 1.3)" model. The red line is prediction of the "Transverse 
Enhancement" model. This color and line style convention is 
used in all subequent plots. Top (a): differential QE cross 
sections. Bottom (b): differential QE cross sections. 



2 

max 



4.9 (GeV/c) 



5.4 Results 

Figures 17] and [8] show the QE differential cross section 
(dcr/dQ^as a function of Q 2 for v^, energies of 1.0 and 
3.0 GeV, respectively. The orange dotted line is the pre- 
diction of the "Independent Nucleon (Ma=1.014)" model, 
the blue dashed line is the prediction of the "Larger Ma 
(Ma=1.3)" model, and the solid red line is the prediction 
of the " Transverse Enhancement" model. The top panels 
(a) show differential QE cross sections, and the bottom 
panels (b) show the Vn differential QE cross sections. 

Figures [9] and [10] show the ratio of the predictions 
of the two models to the predictions of the " Independent 
Nucleon (Af J 4=1.014)" model as a function of Q 2 for v^, 



energies of 1.0 GeV, and 3.0 GeV, respectively. The blue 
dashed line is the ratio for the "Larger Ma (M j 4=1.3)" 
model. The red line is the ratio for the "Transverse En- 
hancement" mode (with error bands shown as dotted red 
lines). The top (a) panels shows the ratio for da/dQ 2 for 
Vp. The middle (b) panels shows the ratio for da/dQ 2 for 
Dfj,. The bottom (c) panels shows the ratio of predicted 
ratio of v^/v^ da/dQ 2 cross sections for the two models 
(divided by the v^/v^ ratio predicted by the "Independent 
Nucleon (Af A =1.014)" model). 

For Q 2 < 0.6 (GeV/c) 2 the differential QE cross sec- 
tion for the "Transverse Enhancement" model is close to 
the "Larger Ma (A/a=1.3)" model. The maximum acces- 
sible Q 2 for 1 GeV neutrinos is 1.3 GeV/c) 2 (as shown in 



figure 16 1. Therefore, fits to the neutrino differential QE 
cross sections for an incident energy of 1 GeV (e.g. Mini- 
BooNE) would yield M A « 1.2 GeV. The extracted value 
of Ma depends on the specific model parameters that are 
used for Pauli blocking and the variation of the statistical 
errors in the data with Q 2 . For a neutrino energy of 1 GeV, 
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Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 



Fig. 9. The ratio of the prediction of the two models for Fig. 10. Same as figure |9jfor Up, energies of 3.0 GeV (max- 
the QE differential cross section da/dQ 2 to the prediction imum accessible = 4.9 (GeV/c) 2 ). 

of the "Independent Nucleon (Ma = 1.014)" model as a func- 
tion of Q 2 for Ufi, energies of 1.0 GeV (maximum accessible 
Qmax = 1-3 (GeV/c) 2 ). The blue dashed line is the ratio for 
the "Larger Ma (Ma=1-3)" model. The red line is the ratio 
for the "Transverse Enhancement" model (with error bands 
shown as dotted red lines). Top (a): ratio for differential 
QE cross sections. Middle (b): ratio for differential QE cross 
sections. Bottom (c): The J? M /^ M ratio for the differential QE 
cross sections divided by the corresponding v^/v^ ratio for the 
"Independent Nucleon (Ma=1.014)" model). 
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the total integrated QE cross section predicted by the the 
"Transverse Enhancement" model is is also larger than 
the total QE cross section prediction of the " Independent 
Nucleon (M A =1.014)" model. 

In the high Q 2 region (Q 2 > 1.2 (GeV/c) 2 ), the pre- 
dicted differential QE cross section for the "Transverse 
Enhancement" model is similar to the prediction of the 
"Independent Nucleon (Myt=1.014)" model. The maxi- 
mum accessible Q 2 for 3 GeV neutrinos is 4.9 GeV/c) 2 . In 
order to reduce the sensitivity to modeling of Pauli block- 
ing, experiments at higher energy [15] typically remove the 
lower Q 2 points in fits for Ma- Consequently, fits for the 
neutrino differential QE cross sections measured in high 
energy experiments would yield a value of Ma which is 
smaller than 1.014 GeV because for Q 2 > 0.5 (GeV/c) 2 
the slope of the differential QE cross section in the tran- 
sition region between low and high Q 2 is steeper than for 
Ma = 1.014 GeV. This is consistent with the fact that the 
average Ma extracted from high energy data on nuclear 
targetsQi] is 0.979 ± 0.016. 

Figure [TTJ shows the total QE cross section as func- 
tion of energy. The data points are the measurements 
from MiniBooNE [B] and NOMAD [IS]. The orange dot- 
ted line is the prediction of the "Independent Nucleon 
(M,i=1.014)" model. The blue dashed line is prediction 
of the "Larger M A (M A =1.3)". The red line is the predic- 
tion of the " Transverse Enhancement" model (with error 
bands shown as dotted red lines). The top (a) panel shows 
the total QE cross section. The middle (b) panel shows 
the total QE cross section. The bottom (c) panel shows 
the ratio of and total QE cross sections. 

Figure [12] shows the ratio of the predictions for total 
QE cross section to the predictions of the "Independent 
Nucleon (Ma=1.014)" model as a function energy. The 
blue dashed line is the ratio of the predictions for the 
"Larger Ma (Ma=1.3)" model, and the red line is ratio for 
the "Transverse Enhancement" model (with error bands 
shown as dotted red lines). The top (a) panel shows the 
ratio of the predictions for the total QE cross section. 
The middle (b) panel shows the ratio of the predictions 
for the total QE cross section. The bottom (c) panel 
shows the predicted v^fv^ cross section ratio divided by 
the predicted v^jv^ ratio for the "Independent Nucleon 
(Myt=1.014)" model. The data points are measurements 
from MiniBooNE [6] and NOMAD [18]. 

As shown in Fig. 11 (a) (top), and Fig. [12] (a) (top) , at 
low energies the " Transverse Enhancement^ model (red 
line) predicts QE cross sections at a level similar to the 
"Larger M A (M j4 =1.3)" model (blue dashed line). Both 
the "Larger M A (M A =1.3)" model and the "Transverse 
Enhancement" model predictions are in agreement with 
the MiniBooNE QE cross sections. However, at higher 
Vfj, energies the " Transverse Enhancement" model predicts 
QE cross sections which are lower than the prediction of 
the "Larger M A (M,a=1.3)". The lower QE cross sec- 
tions at high energy are consistent with the NOMAD mea- 
surements (within experimental errors). 

Similarly, for scattering the "Transverse Enhance- 
ment" model predicts total QE cross section which are 



lower than the predictions of the "Larger Ma (Ma=1.3)" 
model as shown in Fig. [TTJ (b) (middle) and Fig. [12] (b) 
(middle). The lower QE cross sections are consistent 
with the NOMAD measurements (within experimental er- 



6 Conclusion 

We parametrize the enhancement in the transverse QE 
cross section observed in QE electron scattering on nu- 
clear targets as a correction to the magnetic form fac- 
tors of bound nucleons. Within models of MEC, MEC 
processes contribute only to the transverse QE response 
function and do not enhance the longitudinal and axial 
response functions. We find that the QE cross sections 
for v^, QE scattering predicted by the "Transverse En- 
hancement" model agree with the MiniBooNE low energy 
neutrino QE cross sections, and are also consistent with 
QE cross sections measured by NOMAD at higher ener- 
gies. 

The simple two parameter parametrization of Q 2 de- 
pendence of the transverse enhancement as a correction to 
the proton and neutron magnetic form factors can easily 
be incorporated into existing Monte Carlo generators [4]. 



At present, 



experiments use the "Large M4" 



model to predict W 1 



Qelastic 



and W? 



elastic 



for 



neutrino QE scattering on nuclear targets. A large increase 
in Ma is contrary to theoretical expectations |15l[TB] . 

The differential and total QE cross sections predicted 
in the "Larger Ma (M^=1.3)" model are similar to the 
predictions of "Transverse Enhancement" model only at 
low v ^ energies. 

The assumption made in the "Transverse Enhance- 
ment" model is that the enhancements in the transverse 
response functions in v^Jv^ scattering are the same as 
measured in electron scattering, and that there is no ad- 
ditional enhancement in the longitudinal or axial response 
functions (as expected in MEC models [23]). Since we only 
uses parameters from electron scattering data, our anal- 
ysis is purely phenomenological, and does not rely on a 
specific MEC model. Because in electron scattering the 
transverse enhancement is only significant at low values 
of Q 2 , the contribution of "Transverse Enhancement" to 
the total neutrino QE cross section is energy dependent, 
thus resolving much of the apparent discrepancy between 
the low energy and high energy neutrino QE cross sections 
on nuclear targets. 

In an earlier publication, Martini, Ericson, Chanfray, 
and MarteaujH] calculated the contribution of meson ex- 
change currents to the differential and total QE cross sec- 
tions for Pfi/fn energies less than 1.2 GeV. In the com- 
parison with our model, we use the Martini et al predic- 
tions with the random phase approximation ( " QE+np-nh 
RPA"). For the range 0.5 < E < 1.2, the predictions of 
Martini et al are similar to the predictions of the " Trans- 
verse Enhancement" model as shown in Figures |13| and 
14 For E < 0.5, the predictions of Martini et al are lower 
than the predictions of the "Transverse Enhancement" 
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Fig. 11. The total QE cross section as function of energy. 
The data points are measurements of MiniBooNE[6 (gray 
stars) and NOMAD 18 j (purple circles). The predictions for 
the "Independent Nucleon (MA=1.024)" model, "Larger M A 
(Ma = 1.3) model", and "Transverse Enhancement model" are 
shown. Top (a): z/ M total QE cross section. Middle (b): to- 
tal QE cross section. Bottom (c): QE v^jv^ total cross section 
ratio.. 



Fig. 12. The ratios of predicted and measured total QE 
cross section to the predictions of the "Independent Nucleon 
(Ma = 1.014)" model as a function energy. The ratios for the 
predictions of the "Larger Ma (Ma=1.3) model" and "Trans- 
verse Enhancement model" are shown. The data points are 
the ratios for the measurements of MiniBooNE 6 (gray stars) 
and NOMAD [TS] (purple circles). Top (a): The ratio for z/ M to- 
tal QE cross sections. Middle (b): The ratio for P M QE cross 
sections. Bottom (c): The Pp/v^ total QE cross section ratio 
divided by the corresponding ratio for the "Independent Nu- 
cleon (Ma=1.014)" model). 
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Fig. 13. Comparison to the QE cross section predicted by the 
"QE+np-nh RPA" MEC model of Martini et al. [24] (Predic- 
tions for this model have only been published for neutrino ener- 
gies less than 1.2 GeV). The predictions for the "Independent 
Nucleon (MA=1.024)" model, "Larger M A (ATa=1-3) model", 
and "Transverse Enhancement model" are shown. The grey 
squares are the predictions of the MEC model of Martini et 
al. [24] The data points are measurements from MiniBooNE 6 
(grey stars). Top (a): total QE cross section. Middle (b): i> M 
total QE cross section. 
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Fig. 14. Comparison to the QE cross section predicted by the 
"QE+np-nh RPA" MEC model of Martini et a/.[24] (Predic- 
tions for this model have only been published for neutrino ener- 
gies less than 1.2 GeV). Shown are the ratios for the predictions 
of the "Larger M A (M A =l-3) model" and "Transverse En- 
hancement model" to the "Independent Nucleon (Ma=1.014)" 
QE cross section as a function energy. The grey squares are 
the ratios for the predictions of the MEC model of Martini et 
al. |24| The data points are the ratios for the measurements of 
MiniBooNE[6] (grey stars). Top (a): The ratio for total QE 
cross sections. Bottom (b): The ratio for QE cross sections. 



model. However, for such low energies, the predictions are 
sensitive to differences in the modeling of Pauli blocking 
in the two models. The predictions of the Martini et al 
model for v^jv^ scattering for energies greater than 1.2 
GeV have not yet been published. 

Figure [15] shows a comparison of the various model for 
a larger energy range (0.1 to 100 GeV). The energy depen- 
dence for the predictions of the transverse enhancement 
model originates from the energy dependence of the max- 
imum accessible Q 2 (Q^ al ) for QE scattering, as shown 



in Fig 16 The lower energies have lower Q max where the 



; max 

transverse enhancement is large, while higher energies have 
a higher Q 2 nax , where the transverse enhancement is small. 
The differential cross section at high energy is almost in- 



dependent of energy, as shown in Fig. 17 18 19 and 20 



7 Testing the model with neutrino data 

The MINERvA high statistics neutrino experiment [M] at 
Fcrmilab is currently taking data with a fully active scin- 
tillator target calorimeter in the NUMI beam (with both 
neutrinos and antineutrinos) . The QE differential cross 
sections would be measured as a function of Q 2 at a vari- 
ety of neutrino energies within one single experiment, and 
compared to the predictions of various models. 



8 Appendix: v^v^ nucleon/nucleus scattering 

At a fixed value of the final state invariant mass W, the 
differential cross section for scattering at incident 
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Fig. 15. Comparison of predictions for the u^, total QE 
cross section section at high energies for the "Independent Nu- 
cleon (MA=1.024)" model, the "Larger M A (M A =l-3) model", 
the "Transverse Enhancement model", and the "QE+np-nh 
RPA" MEC model of Martini et ai. [24] (Predictions for this 
model have only been published for neutrino energies less than 
1.2 GeV). The data points are the ratios for the measurements 
of MiniBooNE [6] (gray stars) and NOMAD [18] (purple circles) 



energy E is givenJSD] by: 



d(7 G 2 2 „ ?f 1 u , r 21 



+w 2 + w 2 
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2ME 4E ME 



m 2 (Q 2 



AE 2 



m 2 

ME m » 



Here, f£ cos 2 9 C = 80 x lfr 40 cm 2 /GeV 2 . The final state 
muon mass places the following kinematic limits[31] on 
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Fig. 16. The maximum accessible Q 2 for QE events as a func- 
tion of neutrino energy. 

x = Q 2 /2Mv and y = v/E: 

™2 



< X < 1 



2M(£„ - to m ) 
a — b < y < a -\- b , 



(13) 
(14) 



where the quantities a and 6 are 



b = 



1 — 771. 



1 - 



1 



1 



2M£„a; 2E 2 
2 \ 2 2 i ^ 
E 2 , 



/{2 + Mx/E v ) 



2ME„x 



j(2 + Mx/E v ) 



Or alternatively, for a fixed energy and Q 2 , there is a 
maximum value of W which is given bv|32j: 



Wl(Q 2 



1 2 2 

-SO 



4- 2 
777,, mf, 



/ „ m„ 




+ 2 TO m< 



/[a_(Q 2 + m 2 )], 



where s = 2ME+M 2 , a± = l±M 2 /s. For QE scattering, 
this corresponds to a minimum and maximum accessible 
Q 2 for a given neutrino energy. The maximum accessible 
Q 2 {Qmax) for QE events as a function of neutrino energy 
is shown in Fig. 16 



8.1 Quasielastic v^.v^ scattering 

A theoretical framework for quasi-elastic t/ M )-Nucleon 
Scattering has been given by Llewellyn Smith [53]. Here, 
we use the notation of Llewellyn Smith (except that F v 
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Fig. 17. Same as figurej^jfor energies of 10.0 GeV. 



in our notation is equal to £,i s F v in Llewellyn Smith's no- 
tation, where £; s = — 1 — Hn))- In addition, we use Q 2 
while Llewellyn Smith uses q 2 where 



2 2 



4E E' sin 2 - 



-Q 



by 



The hadronic current for QE , v„ scattering is given 



< P(P2)\Jx\ n (Pl) >- 



7a-W) + 



ioj^fJ^ (g 



2Af 



2^ , qxibFpjq 2 ) 



M 



where q = k v — k^, and M = (m p + m n )/2. Here, \x p 
and ^„ are the proton and neutron magnetic moments. 
We assume that there are no second class currents, so the 
scalar form factor T v and the tensor form factor T\ need 
not be included. Using the above current, the QE cross 
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Fig. 18. Same as figure [7] for energies of 25.0 GeV. 



section is 



da v > v 
dQ 2 



M 2 G 2 F cos 2 e c 



8ttE 2 



A(Q 2 )T 
where s — u = AME : 

A(Q 2 ) -- 



(s-u)B(Q 2 ) 
M 2 



C(Q 2 )(a-u) 2 
M 1 



Q 2 



M 2 

+ r (1 - r) \T 2 



f ml 



(l + r)\F A \ 2 -(l-r) \F?\ 2 
(l^f + ^T) 



M 2 AM 2 
{F A + 2F P f - 4(1 + t)T 2 p 



(15) 
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Fig. 19. Same as figure [9] for v^^v^ energies of 10 GeV. Fig. 20. Same as figure [9] for Vy^Vp energies of 25 GeV. 



B{Q 2 ) = ^tT a {TX + ?X) = ItFaGm, (16) Where r = Q 2 /4M 2 . The form factors (Q 2 ) and F^(Q 

are given by: 



C(Q 2 



1^ 



V|2 







(|^| 2 + |F y (Q 



2\|2\ 



(17) 



^W 2 ) + ^aW) 



AM 2 
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Gm(Q 2 ) - G E (Q 2 ) 
+ AM 2 



From conserved vector current (CVC) Q E (Q 2 ) and (jj^ (Q 2 ) 
are related to the electron scattering form factors G E (Q 2 ), 
G n E (Q 2 ), G P M (Q 2 ), and G n M (Q 2 ): 

Ge(Q 2 ) = G p e (Q 2 )-G e (Q 2 ), 

gUq^^gi^-giaq 2 ). 



We also define 



[G v E {Q 2 )] 2 + r[G v M {Q 2 )f 

l + T 



The axial form factor T a can be approximated by the 
dipole form 

2^ _ 9 A 



1 



Q 

M\ 



Where g A = -1.267. 

The pseudoscalar form factor Tp is related to Fa by 
PCAC and is given by: 



2M 2 T A (g 2 ) 



Q 2 



O 2 

Wu— axial c/ \ /-i 



0l^(Q 2 )| 2 



w. 



v— vector 
2 — Qelastic 



w. 



is— axial 
2 — Qelastic 



5{v- 
5{v 



Q 2 ^ 

2M 
2M 



)\Fv(Q 2 t 



)\?a{Q 2 )\ 2 



^^3—Qelastic 



8{u 



2M 



)\2g v M (Q 2 )F A {Q 2 



w 



v— axial 
4—Qelasti 



Q\ i 

) X -X 

2M' 4 



Q 2 



F 2 a(Q 2 ) + 4)\MQ 2 )\ 2 ~ (^a(Q 2 ) + 2F P {Q 2 )) 2 



vector 
5 — Qelastic 



5{v 



v— axial 



i 



5-Qelasttc = 5 ( V ~ QM^^^ 



The vector part of W4 and W5 are related to the vector 
part of W2 and W\ by the following expressions |32|: 



yyvector 
yyvector 



yyvector 



2, ,2 



M 2 V 



Q 4 

.Mv 



.M 2 

o 2- 



Note that: 

a v T ector oc r|^(Q 2 )| 2 ; a a T xlal oc (1 + r)\T A (Q 2 )\ 2 



cx {Ge(Q 2 )) 2 - 



^axial r< 



In the expression for the QE cross section, Fp(q 2 ) is mul- 
tiplied by (ra^/M) 2 . Therefore, in j/„, za, interactions, this 
effect is very small except at very low energy, below 0.2 GeV. 
In the dipole approximation, 

GlAQ 2 )- 4.706 G V D (Q 2 ). 

In our analysis we apply BBBA2QQ7 2 5 corrections [TU] to 
the dipole parametrization of the electromagnetic form 
factors as described in reference [TO]. 

By comparing equations [12] and [15] and using the fol- 
lowing expressions: 

^(Q 2 ) + F 2 V (Q 2 ) = GL(Q 2 ), 

\Tr(Q 2 )\ 2 + T\^(Q 2 )\ 2 = \T v (Q 2 )\ 2 , 

we obtain the following relationships between the struc- 
ture functions and form factors for , QE scattering 
on free nucleons: 



Therefore, for QE and scattering only Qj f con- 
tributes to the vector part of the transverse virtual boson 
absorption cross section. 
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